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Abstract 

Vortex flow in driven type II superconductors shows strong memory and 
history dependent effects. Here, we study a schematic microscopic model of 
driven vortices to propose a scenario for a broad set of these kind of phenom- 
ena ranging from "rejuvenation" and "stiffening" of the system response, to 
"memory" and "irreversibility" in I-V characteristics. 
PACS numbers: 05.50.+q 75.10.-b 47.32.Cc 74.60.Jg 

An important discovery in the dynamics of vortices driven by an external current in type 
II superconductors is the presence of strong "memory" and history dependent effects in 
vortex flow (see for instance and ref.s therein). These off- equilibrium phenomena, par- 
ticularly important at low temperatures, crucially affect the overall system behaviour and, 
thus, have also important consequences for technological applications. Here we present a 
microscopic model to explain their origin and their connections with other "aging" phenom- 
ena of vortex matter such as the rearrangement of vortex domains during creep. In relation 
to recent experimental results we discuss in particular the nature of "memory" effects 
observed in the response of the system to an external drive, i.e., the I-V characteristic. Our 
model explains the peculiar form of "memory" observed in vortex flow at finite T, which 
we call an "imperfect memory". In the same framework one can explain other "anomalous" 
properties such as those of time dependent critical currents. The essential step is to identify 
the relevant time scales in the dynamics. Finally, we stress the similarities and the important 
differences with other "glassy" systems, such as supercooled liquids or disordered magnets 
(see H[10| and ref.s therein). 
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We consider a schematic coarse grained model introduced to describe anomalous relax- 
ation of vortices in superconductors, namely a system of repulsive particles wandering in 
a pinning landscape in presence of an external drive. The model describes several phe- 
nomena of vortex physics, ranging from a reentrant phase diagram in the (B, T) plane, to 
the "anomalous second peak" in magnetisation loops (the "fishtail" effects), glassiness and 
"aging" of slow magnetic relaxation, the "anomalous creep" at very low temperatures, and 
many others ||. 

A system of straight parallel vortex lines, corresponding to a magnetic field B along 
the z-axis, interacts via a potential []J: U(r) = -^frz [K (r / \') — K (r /£')], K being the 
MacDonald function, £ and A the correlation and penetration lengths (£' = c£/\^2, A' = cA, 
c = (1 — B/B^)" 1 ' 2 ). The typical high vortex densities and long A imply that the vortex 
system is strongly interacting. To make it theoretically more tractable, as proposed in 
0]8[], one can coarse grain in the x?/-plane by introducing a square grid of lattice spacing 
I of the order of the London length, A. The number of vortices on the z-th coarse grained 
cell, rii, is an integer number smaller than = B&l 2 /^ || (B C 2 is the upper critical 
field and 0o = hc/2e is the flux quantum). The coarse grained interaction Hamiltonian 
is thus ||: Ti — \ Y^ij n i^ij n j — \ Hi^a n i ~ J2iA P rii. The first two terms describe the 
repulsion between the vortices and their self energy, and the last the interaction with a 
random pinning background. For sake of simplicity, we consider the simplest version of Ti: 
we choose An = A = 1; A^ = Ax < A if i and j are nearest neighbours; A^ = otherwise; 
the random pinning is delta-distributed P(A P ) = (1 — p)S(A p ) +p5(A p — Aq). 

In analogy with computer investigation of dynamical processes in fluids, the time evolu- 
tion of the model is simulated by a Monte Carlo Kawasaki dynamics |TTJ on a square lattice 
of size L [12] at a temperature T. The system is periodic in the ^/-direction. The two edges 
parallel to the y-axis are in contact with a vortex reservoir, i.e., an external magnetic field, 
of density N ext . Particles can enter and leave the system only through the reservoir. The 
above model, called ROM, is described in full details in ||. 

The system is prepared by zero field cooling and then increasing N ext at constant rate 



3 



up to the working value (here, N ext = 10). Then we monitor the system relaxation after 
applying a drive, I (due to the Lorentz force), in the ^/-direction. As in similar driven 
lattice gases [[HJ , the effect of the drive is simulated by introducing a bias in the Metropolis 
coupling of the system to the thermal bath: a particle can jump to a neighbouring site with 
a probability min{l, exp[— (A7Y — eI)/T]}. Here, A7i is the change in 7i after the jump and 
e = ±1 for a particle trying to hop along or opposite to the direction of the drive and e = 
if orthogonal jumps occur. A drive / generates a voltage V []14 |: 



V(t) = (v a (t)) (1) 

where v a {t) = ^ v(t')dt' is an average vortex "velocity" in a small interval around 
the time t. We consider such an average to improve the statistics on V(t) and choose At 
accordingly. Here, t is the Monte Carlo time (measured in units of single attempted update 
per degree of freedom), v(t) = jJ2i v i{t) is the instantaneous flow "velocity", Vi(t) = ±1 if 
the vortex i at time t moves along or opposite to the direction of the drive / and v j = 
otherwise. The data presented below are averaged over up to 3072 realizations of the pinning 
background. 

To characterise the properties of history dependent effects, we consider a striking man- 
ifestation of "memory" observed in experiments where the drive is cyclically changed ||. 
A drive / is applied to the system and, after a time t±, abruptly changed to a new value 
Ii, finally, after waiting a time t 2 , the previous / is restored and the system evolves for a 
further t 3 (see lower inset of Fig.|l]). The measured V(t) is shown in the main panel of Fig.|l| 
for T = 0.1. A first observation is that after the switch to I\ the system seems to abruptly 
reinitiate its relaxation approximately as if it has always been at I\ (see for example the 
dashed curve in Fig.|TJ) , a phenomenon known as "rejuvenation" in thermal cycling of spin- 
glasses and other glassy systems |]10| . The more surprising fact is, however, that for I\ small 



enough (say l\ <C /*, /* to be quantitatively defined below) when the value / of the drive 
is restored the voltage relaxation seems to restart from where it was at t\, i.e., where it 
stopped before the application of l\ (see Fig.|l]). Actually, if one "cuts" the evolution dur- 
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ing t 2 and "glues" together those during t\ and t 3 , an almost perfect matching is observed 
(see upper inset of Fig.[l]). What is happening during t 2 is that the system is trapped in 
some metastable states, but not completely frozen as shown by a small magnetic, as well as 
voltage, relaxation. These non trivial "memory" effects are experimentally found in vortex 
matter H and glassy systems PD |. We call them a form of "imperfect memory", because 
they tend to disappear when the time spent at Ii becomes too long or, equivalently (as 
explained below) when, for a given t 2 , I\ becomes too high, as clearly shown in the inset of 
Fig.[l]. These "rejuvenation and memory" effects are very general in off-equilibrium "glassy" 
dynamics and can be found well above the ideal glass transition point, T c ||, whenever 
the observation time scales are shorter than the intrinsic relaxation scale Ty (to be defined 
below) . 

We now turn to some additional aspects concerning the time dependent properties of the 
current- voltage characteristic. As in real experiments on vortex matter || , we let the system 
undergo a current step of hight Jo for a time to before starting to record the I-V by ramping 
J, as sketched in the inset of Fig||. FigPJ shows (for T = 0.1) that the I-V depends on the 
waiting time to- The system response is "aging": the longer to the smaller the response, a 
phenomenon known as "stiffening" in glass formers P,pD|. These effects are manifested in 
a violation of time translation invariance of the two times correlation functions, which have 
dynamical scaling properties || analogous to those of glass formers |J. 

These simulations also reproduce the experimentally found time dependence of the crit- 
ical current ||. Usually, one defines an effective critical current, 1% , as the point where V 
becomes larger than a given threshold (say V t h r — 10 -5 in our case): one then finds that l e J^ 
is to and Iq dependent (like in experiments 1^ is slowly increasing with t , see Fig.||). 

It is interesting to consider another current cycling experiment which outlines the con- 
current presence of irreversibility and memory effects. The I-V is measured by ramping / up 
to some value l max . Then / is ramped back to zero, but at a given value I w the system is let 
to evolve for a long time t w . Finally, / is ramped up again (see inset of Fig.^). The resulting 
irreversible V(I) is shown in FigfJ. For I > I w the decreasing branch of the plot (empty 
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circles) slightly deviates from the increasing one (filled circles), showing the appearance of 
irreversibility. This is even more apparent after t w : for I < I w the two paths are clearly 
different. Interestingly, upon increasing / again (filled triangles), V(I) doesn't match the 
first increasing branch, but the latest, the decreasing one: in this sense there is coexistence 
of memory and irreversibility. Also very interesting is that by repeating the cycle with a 
new I w (squares), the system approximately follows the same branches. This non-reversible 
behaviour is also found in other glassy systems []10|| . However, spin glasses, for instance, 
show the presence of the so called chaos effects [|K|]. The chaos effect is absent in our system 
as it is also in other ordinary glass formers ||10|| . This kind of interplay between irreversibil- 
ity and memory can be checked experimentally in superconductors and thereby assess the 
present scenario. 

The above experiments outline a set of very important "history" dependent phenomena 
("rejuvenation", "imperfect memory", "stiffening" and "irreversibility" of response, etc.). 
A crucial theoretical step to understand them is the identification of the characteristic time 
scales of the driven dynamics. This we now discuss. 

Upon applying to the system a small drive, /, its response, V, relaxes following a pattern 
with two very different parts: at first a rapidly changing non-linear response is seen, later 
followed by a very slow decrease towards stationarity (see V(t) in Fig.^ for T = 1 and 
/ G {1, 2, 3}). For instance, for / = 3 in a time interval At ~ 2 ■ 10 _1 , V leaps from about 
zero to AVi ~ 2 • 10~ 3 , corresponding to a rate r, = AVi/At ~ 10~ 2 . This is to be compared 
with the rate of the subsequent slow relaxation from, say, t = 2-10 -1 to t — 10 4 , r/ ~ — 10~': 
rj and rf differ of 5 orders of magnitude. 

In agreement with experimental findings |p],^,|5|JT5[| , the slow relaxation of V(t) has a 
characteristic double step structure, which asymptotically can be well fitted by stretched 
exponentials fllEf : V(t) oc exp(— t/ryf. The above long time fit defines the characteristic 
asymptotic scale, ry, of relaxation. The exponent (3 and Ty are a function of /, T and N ext 
(see inset Fig.[|): in particular Ty(I) decreases with I and seems to approach a finite plateau 
for I < I*, with I* ~ 0(1). In this sense, the presence of a drive I makes the approach to 



6 



stationarity faster and has an effect similar to an increase in T. A reasonable fit for 7y(J) 
is: r v (I) = (r£ -r£°)/[l + (J/3/*) a ] +r£°, with t v = 2590, r£° = 1720, I* = 0.86 and a = 2. 

The outlined properties of ry clearly explain the history dependent effects in the ex- 
periments previously considered. For instance, the "imperfect memory" , discussed in Fig.[I], 
is caused by the presence of a long, but finite, scale ry in the problem: for a given I\ the 
system seems to be frozen whenever observed on times scales smaller than Ty(ii). Thus, if t 2 
is short enough (t 2 < ry(/i)) the system preserves a strong "memory" of its state at t±. The 
weakening of such a "memory" found for higher currents I\ in Fig.0, is also a consequence 
of the strong decrease of Ty(I) with /. The phenomenon of "rejuvenation" (see Fig.p]) is, in 
turn, a consequence of the presence of the extremely fast first part of relaxation found in 
V(t) upon applying a drive and of the above long term memory. The existence of the slow 
part in the V(t) relaxation also affects the "stiffening" of the response in the I-V of Fig.|2], 
which is due to the non-stationarity of the vortex flow on scales smaller that ry. Actually, 
in Fig.|2|, for a given / the value of V on the different curves corresponds to the system being 
probed at different stages of its non-stationary evolution (this also outlines that the proper 
definition of I c is the asymptotic one given later on). Finally, in brief, the fact that Ty(I) is 
smaller at high currents, J, and larger at small / (and T), is responsible for the surprisingly 
concomitant effects of irreversibility and memory of Fig|| 

The origin of these time dependent properties of the driven flow, and in turn those of 
I-V's, traces back to the concurrent vortex creep and reorganisation of vortex domains. In 
fact, both with or without an external drive, the system evolves in presence of a Bean like 
profile (see Fig||) which in turn relaxes following stretched exponentials || with a charac- 
teristic time scale tm- An important discovery is that tm{I) and Ty(I) are approximately 
proportional, as shown in Fig.|5|. This outlines that the strong non-linear, non- stationary 
voltage relaxation is structurally related to the reorganisation of vortices during the creep 
(a fact confirmed by recent experiments ||). By assuming that the system is composed of 
evolving vortex domains characterised by drive dependent (free) energy barriers which grow 
with the sizes of the domains, a hierarchy of time scales naturally arises. The simultane- 
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ous presence of both long time scales ("aging modes") and short time scales is a simple 
mechanism able to describe the above phenomena (see ||10|| ). 

Finally, let us notice that the existence at T = 1 of a finite ry limit for I — > (see 
Fig.^) implies that the flux flow can be activated for any finite /. This shows that the 
system's asymptotic (i.e., for t — > oo) critical current, I c , is zero. It is possible that the 
asymptotic value of I c can be non-zero for real three dimensional systems if vortex line 
cutting is negligible. The above also points out that a true divergence in Ty(J), i.e., a true 
asymptotic finite I c , can be found only for T < T c , T c being the ideal glass transition point 
fHH, where by definition r M {T)\ I=0 diverges |§. 

The few key "experiments" we discussed in the context of our schematic microscopic 
model depict a new clear scenario for "aging" and "memory" effects in driven vortex matter. 
In the present framework a broad range of phenomena can be understood (from response 
"stiffening" and "rejuvenation", to time and history dependent critical currents, )■ We 
stress that most of these features can also be found well above the ideal glass transition 
point T c . We have shown how creep and response functions in driven media are related. 
Theoretically intriguing is the observation of similar behaviours found in vortex matter and 
others glass formers such as random magnets and supercooled liquids. 
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FIG. 1. At T = 0.1 for a drive 1=1, the voltage, V(t), is plotted as a function of time. 
As shown in the lower inset, after a time lag t\, the drive is abruptly changed to I\ for a time 
ti and finally it is set back to its previous value. When / is switched to I\ the system seems 
to "rejuvenate": it suddenly restarts its relaxation along the path it would have had if / = I\ 
at all times (consider the continuous and dashed bold curves, corresponding to I = I\ = 1 and 
/ = I\ = 0.8, plotted for comparison). By restoring / after t2, the system shows a strong form of 
"memory": if ti and I\ are small enough (see text) the relaxation of V(t) restarts where it was at 
ti. However, if ti and I\ are too large, this is not the shown in the upper inset. In this 

sense, the above is an "imperfect memory". 




11 



> 



0.0014 
0.0012 
0.001 
0.0008 
0.0006 
0.0004 
0.0002 




50 100 150 200 
t 



■ t =10' 

A t =10 

♦ t =10 

• t n =10 7 



!•••• 







0.0 0.5 1.0 1.5 2.0 2.5 3.0 

I 



FIG. 2. The I-V obtained at T = 0.1 by ramping / after keeping the system in presence of a 
drive /o = 1 fo r a time to as shown in the inset. The response, V, is "aging" (i.e., depends on to) 
and, more specifically, stiffening: it is smaller the longer to- 
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FIG. 3. The I-V is measured at T = 0.1 during cycles of / (see also the inset): I is at first 
increased up to I ma x (filled circles) ; along the descending branch of the cycle (empty circles) , when 
I = I w (in the main panel the I^'s, for two cycles, are located by the arrows) the drive is kept fixed 
for a time t w = 10 4 and then the cycle restarted; finally, / is ramped up again (filled triangles). 
For I > I w , the first increasing ramp and the decreasing one (resp. filled and empty circles) do 
not completely match, showing irreversibility in the I-V. After waiting t w at I w , a much larger 
separation is seen. However, by raising / again (filled triangles) a strong memory is observed: 
the system doesn't follow the first branch (filled circles), but the decreasing one (empty circles). 
Furthermore, in a cycle with a lower I w (squares), the same branches are found. 



13 




-2-10 1 2 3 4 5 

10 10 10 10 10 10 10 10 

t 



FIG. 4. The time evolution of the response function, V(t), for the shown values of the drive / 
(at T = 1 and N ex t = 10). In the asymptotic regime V(t) is well fitted with: V(t) oc exp[— (t/ryf]- 
Inset The characteristic scale of relaxation, 7y(J), as a function of /. For I — * 0, Ty(J) seems to 
saturate to a finite value which implies I c = 0. 
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FIG. 5. Left Vortex driven flow and creep are strictly related. Actually, the characteristic 
time scales of voltage relaxation, Ty(I), and of magnetic creep, tm{I) (shown for T = 1), are 
approximately proportional. Right The Bean field profile, N(x), present in the system during 
relaxation is shown for the plotted values of I and t. 



